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This study presents a control scheme of the electronic interface of a grid 
connected Variable Speed Wind Energy Generation System (VS-WEGS) 
based on a Doubly Fed Induction Generator (DFIG). The efficiency of the 
wind energy is represented to according to the control strategy applied. Thus, 
in the case of unbalanced grid voltage, the negative sequence voltage causes 
additional strong oscillation at twice the grid frequency in the stator 
instantaneous active and reactive powers. The objective of this work is to 
present an enhanced MPPT controller uses backstepping approach 
implemented in both dq + and dq“ reference frames rotating to keep a safe 
operation of DFIG during unbalanced grid voltage associated with regulating 
rotor flux, in order to estimate rotor flux, a nonlinear observer based on 
sliding mode is proposed in this work. Note that the conventional controllers 
(PI, PI-R ...) are not provided satisfactory performance during unbalanced 
voltage dips. The validation of results has been performed through simulation 
studies on a 4 kW DFIG using Matlab/Simulink 
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1. INTRODUCTION 

Wind power system can be exploited at constant speed or variable speed operations using power 
electronic converters. In this work, the wind turbine will be exploited at the maximum power operating point 
(MPPT) for various wind speeds by optimally controlling the speed of the shaft [1]. In contrast, regulating 
the rotor flux need the state vector of the controlled system.Thus, a good knowledge of the rotor flux is a 
necessity [2]-[3]-[4].Hence; the use of a physical sensor introduces a lot of problem in maintenance, quality 
of service, precision and cost. To overcome this difficulty, preceding works present the observer for flux 
estimation, then the speed is estimated from flux [5] - [6] .Indeed; different models have been proposed for 
estimating rotor flux. One distinguishes for example statistical models (e.g. [7]). The proposed estimators 
generally show a complex processing mode and, therefore, lead to serious implementation problems in real 
time. For this purpose, the Sliding Mode Observer presented in this article will be based on measuring 
machine currents and voltages only. 

However, the unbalanced network voltages have detrimental effects on the operation of the DFIG 
when the stator is directly connected to the grid as in our case study [8]. In this sense, and due to the strong 
penetration of DFIG into wind turbine structures, much research in recent years has sought to improve robust 
and accurate control under various network disturbances [9]-[10]-[ll]. With the increasing use of wind 
power, it has become unacceptable to accept disconnection of the wind turbine with each occurrence of a 
fault on the electrical grid. This disconnection could result in a large number of undistributed energy, which 
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favors network instability. However, the grid voltage has been assumed symmetrical for many works. While 
in practice, asymmetric faults are often present. For a DFIG generator, if the unbalanced voltage has not been 
taken into account in the control law, with a shallow depth, high pulses can be obtained at twice the grid 
frequency appearing in the torque electromagnetic, active and reactive powers of the stator, which implies the 
separation of the wind turbine system of the grid [9]-[12]-[13]. In [14]-[15] the PI controllers are used to 
control the DFIG under unbalanced condition. Some authors use PI-R controller where the positive 
synchronous reference frame can directly regulate both positive and negative sequence components without 
involving sequential decomposition [16]. Moreover, some strategies such as DPC based on power definition 
connected to unbalanced grid voltage are proposed with reduced double grid frequency oscillations of 
electromagnetic torque [17]-[18], but despite the simplicity of these controllers, they are not enough at the 
point of precision. 

In this work, the performance of the wind integration under unbalanced grid has been enhanced by 
the elimination of oscillations indueced by the asymmetric voltage, based on rotor flux sliding mode 
observer. The paper is organized as follows: the wind turbine modeling under unbalanced grid voltage is 
presented in section 2. Section 3 is devoted to design of the rotor flux sliding mode observer and controllers’ 
laws. The control performance is shown by simulation in Section 4, a conclusion and reference list end the 
paper 


2. WIND TURBINE MODEL UNDER UNBALANCED GRID 
2.1. Turbine model 

The turbine model is inspired from [1]. The available aerodynamic wind power is given by: 

Paero=JpJlR 2 Cp(A,3)v^ (1) 

Where (3 is a blade pitch angle, in this paper, (3 is held constant (C P (A, (3) = C P (A)), C P (A) is the 
coefficient of performance, v w is the wind velocity, p is the air density, R is the radius of the wind 
turbine rotor and the tip speed ratio A is given by 


2.2. Dynamic DFIG model with iron losses 

The equivalent circuit of a DFIG is inspired from the proposed model in [19]-[20] .When the iron loss 
resistance R m is connected in parallel with the stator inductance (L sa and M). The choice of this model 
is justified mainly by the simplicity of calculations, because we will have a similar circuit to the 
conventional model (without iron losses) using thevenin’s theorem. The equivalent circuit of the DFIG 
after arrangement is shown in Figure 1. In addition, the R m resistance is a variable parameter.The 
Thevenin’s equivalents for stator resistance and voltages/currents are calculated in the stationary, a- (3 
reference frame as follows [19]-[20]-[21]. 
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( 3 ) 




Figure 1. DFIG equivalent circuit in the dq + reference frame: (a) direct axis, (b) quadrature axis [19] 
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The model considered for DFIG can be written as follows: 

f x = f(x,u,w) = [fiix.u, w) f 2 (x,u, w) ... f s (x,u, w)Y with \x = [i+ Td if Tq cpf d (pf q n] T 
fy = - Vh(x) h 2 (x) h 3 (x)] T \u = [ vf Td vf Tq ] T and w = 

[ Vrd Y 




■< 


fl(x,U, w) = — r|ij T d + « s isTq + aP^Prf + P (J J<Pr + q + ^ V s + Td - pv r + d 

/ 2 (x, u, w) = -co s if Td - r]is Tq - Pu(p? d + ap<p? q + 2-V s + Tq - (3V r + q 

f 3 (x,u, w) = aMi+ Td - a(pr d + gM s <P?q + K + d 
Uix.u.w ) = aMi^ Tq - go) s (Prd ~ a<p? q + V r + q 
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( 4 ) 


H : DFIG rotor speed; J : total inertia constant; f : viscous friction coefficient; (L s , L r , M): 
Stator, rotor and mutual cyclic inductance; (L sa , L rcj ): stator ,rotor cyclic leakage inductance ; (R s , R r ): 
stator and rotor resistances ; (cp s , cp r ): stator , rotor flux components ; Superscripts (+, -): positive, negative 
reference frame ; Subscripts (+, -) : positive, negative sequence component. 

Figure 2 presents the spatial relationships between the stationary (ap)s reference frame, the rotor 
(aP)r rotating at the speed of oo r , and the dq + and dq _ reference frames rotating at the angular speed of co s 
and - 0 ) s , respectively. 

fd + q = ; /- = f {ap)s e^st ;/+ = f {af 0 r e-J°>s Up+ t . f - = f ^ e -jco sliv .t 

Where: f represents the voltage (current and flux), oo slip+ = oo s — go and oo slip _ = — oo s — oo 

According to (4), the stator current can be calculated as: 


■ + _ ^sdq M*rdq 

'sTdq - ~ 


From (4) and (5), the rotor flux can be written 
Mc Psdq 


<Prdq 


- + L r CTi^ 


rdq 


Substituting (6) into (4) yields the rotor voltage in the dq + reference frame as 

+ d ■+ . M d^sdq , T ^rdq , ■ /^^sdq . » ■+ \ 

Vrdq = RrCdq + “^ + 4®“** + jm slip+ * + L r CTl rdq J 


( 5 ) 


( 6 ) 


( 7 ) 



Figure 2. Relationships between stationary (aP)s reference frame, the rotor (aP)r reference frame 
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2.3. Instantaneous active and reactive power flow under unbalanced grid 

By neglecting the stator resistance voltage drop and considering (4), the stator voltage can be 
represented in the positive dq + reference frame as [19] 


^sTdq 


jo> s (cp+ q+ - cp s “ dq _e H2 “ st ) 


( 8 ) 


According to (8), the stator flux in the positive dq + reference frame can be expressed evidently as: 

- i 2a) s t (9) 


< Ps + dq=“( V W“ V ^d q - 


Under unbalanced grid voltage conditions, the instantaneous active and reactive power outputs from 
DFIG stator can be expressed in (10 ) with T s xd q is the conjugate complex of i s j dq : 


, -+ 


P S + jQ s = V sTdq X 'sTdq 


( 10 ) 


Substituting (5), (8) and (9) into (10) and separating the instantaneous active and reactive powers 
into different pulsating components yield [19]: 


Ps PsO "P Ps—sin2 Sin^2u) s t) T P s — CO s2 COs(2(U s t) 
Qs = Qso + Qs-sin2 Sin(2w s t) + Q s _ cos2 cos(2w s t) 
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(11.a) 
(11 .b) 


( 12 ) 


3. ROTOR FLUX SLIDING MODE OBSERVER AND MPPT CONTROL STRATEGY FOR 
DFIG 

The aim of the control is to optimize the extraction of aerodynamic power (MPPT) under an 
unbalanced network voltage. The reference active power is determined as a function of the wind speed. So, 
based on the current model of the DFIG, a sliding mode full state observer is presented to estimate the rotor 
flux [22]-[23]. 

3.1. Rotor flux sliding mode observer 

In this section, we present the rotor flux sliding mode observer for the MPPT control with rotor flux 
regulation. The simplicity and ease of implementation among the advantages of this observer. Figure 3 shows 
the structure of the rotor flux observer. 



Figure 3. Rotor flux sliding mode observer structure 
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Let us rewrite the DFIG model (4) in the following form: 
•i s + T = A 1 1 Jr + + 6u - |3w 


{ 


(13) 


0+ = aMI+ T + B 2 <D+ + w 


(14) 


With ^ = _;];*,=P[_“„ “] '• «2 = [_ 8<u , ®«1^-S;'.'Jr = [* a *,] T and« 

[ Tbd *Prq ] 

Hence, the structure of the observer can be expressed by [22]-[23]: 

ri s + T — + 5 u — ( 3 w + DjU s 

= aMIJr + £ 2 0+ + w + D 0 u s 

The dynamics of the estimation error are expressed by the following equations: 

J ^sT = ^l^sT + — D t U s 

^ $+ = aMl s + T + £ 2 $+ - D 0 u s 


+ _ 

r 


(15) 


With : I<f T — I<f T — I<f T ; 0+ — 0+ — 0+; u s — sign(S 2 )] r ; S — [5^ S 2 ] T — I<f T 

S is the sliding mode surface and (D;; D 0 ) are the matrixes (2x2) that we will determine later. 
Consider the Lyapunov candidate function V s = -S T S , we obtain: 


14 = S T S = S r I s + T = 5 r (^ 1 I s + T + #!$+) - S T D i u s 
in order to satisfy the condition of attractiveness ( S T S < 0) , we must have 
S T (AJ S + T + B 1 ®+) < S T D iUs 


If we put: D t 


Mi 0 
0 M 2 


. Then we obtain the condition below 


hi|S 1 | + q 2 |S 2 |>S T (A 1 I s + T + B 1 O r + ) 

When the sliding mode is reached, the switching surface will verify it: I<f T = I<f T = 0 

Therefore, we obtain: u s = D^B^f 

Introducing (19) in (15), we obtain 

$+ = (B 2 -D 0 D i - 1 B 1 )$+ 

We put: — B 2 +D,pD[ 1 B 1 -P 

Equation (20) becomes: = — P<t>+ 

In order to have exponential convergence, we choose P in the following form: P = 
Where p x and p 2 are positive constants. 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


Pi o 
.0 p 2 


Finally, (21) can be rewritten as follows: D<j, = ( 


Pi o 
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To complete the design of this observer, we must choose the convenable observer parameters. We 
note that the (p x ; p 2 ) parameters determine the dynamics of observer convergence and the (p x ; p 2 ) 
parameters verify the conditions of attractiveness and stability of the observer. 


3.1. Maximum wind power extraction 

In order to capture the maximum power of the incident wind, we must adjust the rotational speed of the 
turbine permanently to the winds. The optimum operation speed of the generator is estimated for A opt = 8,2 

by: 


GV W X. 


W opt - 


opt 


( 22 ) 


Where G is Gearbox ratio of the wind turbine. 


3.2 Backstepping control design under unbalanced grid voltage conditions 

The control objective is twofold: the objective of higher priority is tracking the references trajectory 
(P s _ re f, Qs-ref) where P s _ re f is the optimal stator active power reference according to optimal speed of the 
generator and the second is to keep the rotor flux at its nominal value by controlling the stator reactive power 
given via (27) regardless of the asymmetric voltage. 

Using stator voltage orientation approach, i.e.V s + Td+ = 0, the rotor current references can be 
considered for eliminating the double frequency pulsations of stator output active power (P s - S in 2 = 
P s _cos2-0)- According to (11-a), the references of the positive and negative sequence rotor currents are 
calculated as [19]: 


Wd+_ref 
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~ L s v trq+ 
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2v sTq- 
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(23) 

(24) 

(25) 

(26) 


With: D 1 = vtr q + + v s ? d _ + v j q _ and D 2 = v+ T z q+ ~ (vj d _ + vj q _) 

And P s _ re f = p s0 • Also, equations (6) and (23) indicate respectively that the rotor flux and the stator 
reactive power Q s0 can be controlled by the d+ axis rotor current components (i+d+ re f), using (6)-(23)-(24)- 

(25) and (26), the signal reference (0 re f) of the rotor flux magnitude (Jtyrd + (Prq ^ can be expressed 

according the reactive power reference Q s0 by: 


0re/ MD * QsO + + l rd-) + ^ + [ L r^rq] 


From (27), Q s0 can be easily extracted ( for 0 ref = 0,45 Wb ). 

Using (7), (23), (24), (25) and (26) the errors dynamic equations are given by: 
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(27) 


(28) 

(29) 

(30) 
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■ d / v sTd- 'rd+ +v sT 

4 dt V v sTq+ 


v sTq- 1 rq+ 
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) aL r t Vrc l _ ^rbq- ^slip-^Prd-] 


Consider the Lyapunov candidate function 


V = - e 4 + - e| + - e 3 + - e 4 

2 2 z 2 6 2 


(31) 


(32) 


To attain the objectives of tracking performance, the control signals v+ d+ , v+ q+ , v rd _ and v rq _ can be 
extracted from (28), (29), (30) and (31) satisfying V<0 [19]: 
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(33) 

(34) 

(35) 

(36) 


Where k lf k 2 , k 3 and k 4 are positive constants.Introducing (33), (34), (35) and (36) in (32), the time 
derivative of the candidate Lyapunov function becomes negative definite: V = — k ± e\ — k 2 e 2 — /e 3 e| — 
/c 4 e 4 < 0. This assures the global asymptotic stability of the error system. Based on the proposed control law 
described by (33), (34), (35) and (36), Figure 4 shows the schematic diagram of the control system for a 
DFIG. 


3 <f> supply 



Figure 4. The schematic diagram of the proposed control 


4. SIMULATIONS RESULTS 

In order to demonstrate the performance of the rotor flux sliding mode observer, the simulation 
procedure is designed to consider a constant rotor flux (equal to its nominal value : 0 ref = 0,45 Wb). The 
observer parameters (Pi,p 2 ) and (q^ \i 2 ) are chosen as follows: p x = 200 ,p 2 = 60, = 1000 and |i 2 = 

1000. Figure 5 shows the satisfactory performances of the rotor flux estimation. The estimated rotor flux 
tracks the actual rotor flux with good precision and the estimation error obtained is practically zero. 
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Table I: Turbine parameters 

Table II: DFIG parameters of DFIG 4 kW 


Number of blades 

3 

Number of poles pairs of the DFIG 

2 

Turbine radius, R (m) 

3 

Stator resistance 

1.25 Q 

Gear box ratio 

1:8.5 

Stator leakage inductance 

0.00096 H 

Power 

5 kW 

Rotor resistance 

0.170 

Cut-in wind speed 

4.25 m/s 

Rotor leakage inductance 

0.0018 H 

Rated wind speed 

10.5 m/s 

Mutuelle inductance 

0.0772 H 

Rotor inertia 

0.05 kg m2 

Viscous friction coefficient 
Generator inertia 

0.001 m/s 

0.33 kg m2 



0.6 1 1.6 2 2.6 3 3.6 4 4.6 6 6.6 6 


Figure 5. Rotor flux estimation results. Upper: Actual and estimated rotor flux; lower: estimation 

error 


4.2. Controllers performances 

The stimulation procedure is designed in such a way: during (0s-2s), the conventional Backstepping 
control was established during steady state (under the stable grid). Then at 2s, we started an asymmetric fault 
(a 20% of the nominal grid voltage dip of V sa ) to show the impact of the asymmetric fault on various signals 
of DFIG (with the conventional control). At the end, for (4s-7s) we show the feasibility of the enhanced 
Backstepping control based on the control law described by (33), (34), (35) and (36). The control parameters 
are chosen as follows: k x = 30 , k 2 = 120, k 3 = 115 and k 4 = 500. The satisfactory performances of 
the stator active and reactive power tracking are shown in Figure 6.b and Figure 6.d by limiting the 
mechanical stress induced by asymmetrical fault. Also, note that for imposing a constant rotor flux, the stator 
reactive power reference is imposed according to (27), see Figure 6.c.We see that the error values remains 
practically zero 




(a) Rotor Flux (Wb) 



Figure 6. Simulation responses for 20% of depth of V sa 
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5. CONCLUSION 

In this article, an improved MPPT control is developed using backstepping technique under 
unbalanced grid voltage associated with a rotor flux observer based on the sliding mode approach. The fault 
mode is activated when an asymmetric fault occurs in the grid voltage in order to eliminate the oscillations of 
the signals. Reference rotor currents are calculated, in order to maintain safe operation, as today's grid codes 
require. The proposed rotor flux observer is based on measuring electrical signals (stator currents, rotor and 
stator voltages). The performance and stability of control laws and observer have been shown by simulation 
studies in Matlab/Simulink® environment. 
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